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We consider the most general new physics effective Lagrangian for b — > sl + l~. 

IT) 

O ' We derive the upper limit on the branching ratio for the processes B s — > l + l where 

O ' 

■ I = e,fi, subject to the current experimental bounds on related processes, B — > Kl + l~ 

H : 

and B — ► K*l + l . If the new physics interactions are of vector/axial- vector form, 
the present measured rates for B — ► (if, K*)l + l~ constrain -B(S S — > Z + Z~) to be of 
the same order of magnitude as their respective Standard Model predictions. On the 
other hand, if the new physics interactions are of scalar /pseudo-scalar form, B — > 



> 

\ (K, K*)l + l~ rates do not impose any constraint on B s — > l + l~ and the branching 

fS| ' ratios of these decays can be as large as present experimental upper bounds. If future 

o ' 

experiments measure B(B S — > l + l ) to be > 10 8 then the new physics giving rise 

o 

c-j . to these decays has to be of the scalar/pseudo-scalar form. 

P 



The rare decays of B mesons involving flavour changing neutral interaction (FCNI) b — ► s 
has been a topic of great interest for long. Not only will it subject the standard model (SM) 



^ ■ to accurate tests but will also put strong constraints on several models beyond the SM. In 

H ' 

the SM, FCNI occur only via one or more loops. Thus the rare decays of B mesons will 
provide useful information about the higher-order effects of the SM. Recently, the very high 
statistics experiments at B-factories have measured non-zero values for the branching ratios 



for the FCNI processes B -> (K, K*)l+l~ []J 

Br(B — > Kl + l~) = (4.8^;° ±0.3 ±0.1) x 10" 7 , 

Br(B -> K*l + T) = (11.5±|;| ± 0.8 ± 0.2) x 10" 7 . (1) 

These branching ratios are close to the values predicted by the SM Q]. However, the SM 
predictions for them contain about ~ 15% uncertainty coming from the hadronic form 
factors. Still, it is worth considering what constraints these measurements impose on other 
related processes. 
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The effective Lagrangian for the four fermion process b — > sl + l~ gives rise to the exclusive 
semi-leptonic decays such as B — > Kl + l~ and B — * K*l + l~ and also to purely leptonic 
decays B s — > where I — e, 11. (From here onwards, the symbol I represents either 

e or fi.) Relation between semi-leptonic and purely leptonic B-decays, arising from FCNI 
generated by heavy Z' boson exchange, was briefly considered in J, [5j. The SM predictions 
for the branching ratios for the decays B s — > e + e~ and B s — > n + ji~ are (7.58 ± 3.5) x 10~ 14 
and (3.2 ±1.5) x 1CT 9 respectively jf]|. The large uncertainy in the SM prediction for 
these branching ratios arises due to the 12% uncertainty in the B s decay constant and 10% 
uncertainty in the measurement of V ts . These branching ratios have been calculated in 
various newnhysics models. In models with Z'-mediated FCNI, one has B(B S — > /i + /i~) < 
5.8 x 10~ 8 |l| which is about 20 times larger than the SM prediction. Due to the increased 
precision in the measurement of B{B — > (K, K*)l + l~) 7 this bound can be improved and the 
present calculation attempts to do so. B(B S — > l + l~) are also calculated in multi Higgs 
doublet models. These models are classified into two types. In the first type, there is 
natural flavour conservation (NFC) and there are no FCNI at tree level. In such models, 
there is an additional loop contribution to FCNI, where a charged Higgs boson exchange 
replaces the SM W-exchange. In a two Higgs doublet model with NFC, branching ratio 
for B s —>■ - > !0- 8 1S p „ Ssibl e B. iT*. second t yP e, flavour cflan^ p_ 
do occur at tree level, mediated by flavour changing neutral scalars (FCNS's). In such 
models also a branching ratio of about 10~ 8 for B s — ► can be achieved From the 

experimental side, at present, there exist only upper bounds B(B S — > e + e~) < 5.4 x 10~ 5 
Q and B(B S -> ^^~) < 5.0 x 10" 7 Q. 

In this paper, we consider the most general four fermion effective Lagrangian for b — > sl + l~ 
transition due to new physics. We derive upper bounds on the branching ratios for B s — ► 
e + e~ and B s — > by demanding that the predictions of this new physics Lagrangian for 

B — > K*l + l~ and B — > Kl + l~ should be consistent with the current experimental values. 

The most general effective Lagrangian for b — > sl + l~ transitions due to new physics can 
be written as, 

L eff (b^sl + r) = L VA + L SP + L T (2) 

where, L V a contains vector and axial-vector couplings, L SP contains scalar and psuedo- 
scalar couplings and Lt contains tensor couplings. Lt does not contribute to B s — > l + l~ 
because (0|scr M ^6|i? s (j9 B )) = 0. Hence we will drop it from further consideration. 
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First we will assume that the new physics Lagrangian contains only vector and axial- 
vector couplings. We parametrize it as 



Gp ( Oi 



L VA {b -> sl + l ) = ^= [j^T J s (9v + dAlshu l(g' v + 9a1s)1^ 1 - (3) 

Here the constants g and g' are the effective couplings which charecterise the new physics. 
From the above equation, we get B s — > matrix element to be 



&(- 



M(B S - Z+r) = (g A g A )^ [ j^r (0\s l5ltl b\B s )(l + l- 7767*1 0). (4) 



Only the axial vector parts contribute for both the hadronic and leptonic parts of the matrix 
element. Substituting (0 |s757^&| B s ) = —ifB s PB^, in Eq. (J3J we get 



u{pi)lsv{pi). (5) 



As we are considering only vector and axial vector currents, helicity suppression is still 
operative for the B s — > l + l~ decay amplitude. The calculation of the decay rate gives 

T NP (B S - / + r) = SSl (gA9 A ) 2 m Bs ml (6) 

Thus the decay rate depends upon the value of [gAg A ) 2 . To estimate the value of (<7a5u) 2 ? we 
consider the related decays B — > K*l + l~ and B — > Kl + l~ , which also receive contributions 
from the effective Lagrangian in Eq. (JHJ. In deriving Eq. Q, we dropped terms proportional 
to m 2 /m B , as their contribution is negligible. We will make the same approximation in 
calculating the decay width of semi-leptonic modes also. 

We first consider the process B — > K*l + l~ . Here we will have to calculate the following 
hadronic matrix elements 

(K*(p K *)\s-f^b\B(p B )) = ie^ Xa e v (j> K *){p B +PK*) X {PB-PK*YV{q 2 ) 
(K*{p K ,)\s lhlfJ b\B{p B )) = e^p K ,){m 2 B -m 2 K ,)A l {q 2 )-{e.q){p B +p K ») ll A^q 2 ) (7) 

where q = pi+ +Pi-. In the above equation, a term proportional to is dropped because 
its contribution to the decay rate is proportional to mf/m B . It is assumed that the q 2 
dependence of these form factors is well described by a pole fit: 

V 



V(q 



2\ 



(m B + m K *){l - q 2 /m 2 B ) 



AM 2 ) - A ' 



m B + m K *)(l - q 2 /m 2 B ) 
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The decay rate is 

2 

2 V 192tt 3 / V Ans~ 



TMB KVl') = i f^|) ( ° ] to? + 9 2)/ v x, (8) 



, W, 



where Iva is the integral over the dilepton invariant mass (z = q 2 /nig). Under the assump- 
tion that A\ ~ A 2 , Iva is given by 

PZmax Z r, r, r z max 7 

Iva = g v V 2 / tfe-— Ci(z) + / dz-—C 2 (z), (9) 



where, 



nix* x 2 



Ci(z) =2 1 + -^ $(z) 



3(1- + f-^y fi + ^r 2 (* - ^\ $(.) 



V nig J \lniK* / V 771 5 / V 777^ 

with $(2) = (1 — z) 2 + 4z (m K * /mgf ' . The limits of integration for z are given by z min = 
{2mi/m B ) 2 and z max = (1 — m K */m B ) 2 . From equation (JHJ) we see that, the value of 
{9a9a) 2 can be determined from the measured rate of T(B —> K*l + l~), provided the value 
of gy(g'v + 9a) is known. For this we consider the decay of B — > Kl + l~ . 
The matrix element neccessary in this case is Q] 

(tffc*) |s 7 ^l B{p B )) = ( PB +PK)^ B (q 2 ), (10) 

where again a term proportional to is dropped. The q 2 dependence of the formfactor, 
again, is approximated by a single pole with mass ~ nig, 



The decay rate is given by 



We demand that the maximum value of this decay rate is the measured experimental value, 
(ie.) 

Texp = TiVP- (13) 

With this assumption we calculate the upper bound on the decay rate of B s — ► / + / _ , arising 
due to LvAi given in Eq. (j3J). Using Eqs. (JHJ, (jT2"l) and (fTT^ . we get 

2 / '2 , '2\ Bex V {B Kl + l ) 1n 4 / 1 .s 

QviQv + 94 = 9 x 10 14 

yyvyy ^ 3.45 [/+(0)] 2 1 J 
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and 

, a , _ B Exp (B -» x 10 4 - 1.58V 2 ^(^ 2 + ^ 2 ) 

SUlfiV + - g g^2 • l i0 J 

In our calculation, we take the formfactors to be 

/+(0) = 0.3191°;°!? 
V = 0.457^ 

A, = 0.337^, (16) 

and use experimental values of B — > (K, K*)l + l~ given in |2j. Adding all errors in quadra- 
ture, we get 

fibff + gl) = (1.36l°;2) x io- 2 

si(^ 2 + g'i) = x io- 3 . (17) 

Thus the maximum value (5U5U) 2 can have, is 

{g A g' A f = (6.761J3) x 10~ 3 (18) 
The branching ratio for 5 S — > l + l~ , due to Ly^, to be 

S (S s ^ e + e -) = i.oe x IO" 10 ■ f 2 Bs (g A g' A ) 2 

5(5 S - n+yT) = 4.54 x 10" 6 • fl(g A g' A ) 2 . (19) 

Substituting /g s = 240 ± 30 MeV [12] and the maxmimum value for (g A g' A ) 2 from Eq. JTHJ), 

we get 

B(B S -> e+e") = 4.06±1;H x 10~ 14 

5(5 S -> /i+/0 = 1.74 4 :};! 3 x 10~ 9 (20) 
Therefore the upper bounds on the branching ratios are, 

B(B S -> e + e~) < 6.71 x 10~ 14 

5(5, -> < 2.87 x 10~ 9 (21) 

at la and 

5(5, -> e+O < 1.20 x 10~ 13 

5(5, -> < 5.13 x 10~ 9 (22) 
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at 3<r. 

These rates are close to the SM predictions. The reason for this is quite simple. The 
decay rate for an exclusive semi-leptonic process can be written as 

T = (c.c. ) 2 (/./.) 2 phase space, (23) 

where c.c. is the coupling constant and /./. is the form factor. The measured rates for the 
exlcusive semi-leptonic deays are close to the SM predictions. And we assumed that the 
new physics predictions for these processes are equal to their corresponding experimental 
values. Also, the same set of form factors are used in both SM and new physics calculations. 
Thus the assumption that new physics predictions for semi-leptonic branching ratios are 
equal to their experimental values (which in turn are equal to their SM predictions) implies 
that the couplings of new physics are very close to the couplings of the SM. This is why 
our new physics prediction for the purely leptonic mode is also close to the SM prediction. 
Therefore, new physics, whose effective Lagrangian for b — > sl + l~ consists of only vector and 
axial vector currents, cannot boost up the rate of B s — > l + l~ due to the present experimental 
constraints coming from the decays B — > Kl + l~ and B — > K*l + l~ . 

For the reasons explained above, using a different set of form factors, as for example those 
given in [uj, will not change the upper bound on B s — > l + l~ significantly. In fact, we find 
that the change is less than 10%. 

We can obtain a more stringent upper bound on {gAQA) 2 by the following procedure. We 
equate the new physics contribution for T(B — » (K, K*)l + l~) to the difference between the 
experimental value and the SM contribution. This, in turn, leads to a much more stringent 
upper bound on contribution of Lva to B s — > l + l~ . In fact, at la, this bound is consistent 
with 0. At 3<T we get 

B(B S -»■ e+e~) < 7.89 x KT 14 

B(B S -> fi+fi-) < 3.37 x 1(T 9 , (24) 

which are again comparable to the SM predictions. Comparing these results with the ones 
obtained by previous assumption, we see that there is not much difference in the branching 
ratios. This occurs due to the relatively large errors in both the experimental measurements 
and SM predictions for T(B — » (K, K*)l + l~). Thus we conclude that the presently measured 
values of B — > (K, K*)l + l~ do not allow any large boost in the contribution of Lva to 
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We now consider the new physics effective Lagrangian to consist of scalar/pseudoscalar 
couplings, 



L SP (b sl + l ) = ^| (^j^r^j s (9s + gpis)bl{g' s + g'pisil- 

The matrix element for the decay B s — > l + l~ is given by, 



On substituting, 



a 



V2 \^s 2 w/ 

(0 |s 75 6| B B ) = -i 



g P (0 |s7s&| B s ) g s u(pi)v(pi) + g P u[pi)^v[pi) 



we get, 



M(B s ^l + r) = -ig p ( ^=( 



a 



fB s m 2 Bs 

A-Ksly J m b + m s 



fB s m 2 Bs 
m b + m s 



g' s u(pi)v(pi) + g' P u(pi)~f 5 v{pi 



(25) 



(26) 



(27) 



(28) 



where m b and m s are the masses of bottom and strange quark respectively. Here we take 
the quark masses from Particle Data Group obtained under MS scheme [lJ]. We see that 
in this case there is no helicity supression i.e. the rates for the decays B s — > e + e~ and 
B s — > will be the same provided g' P and g' s for both electrons and muons are the 

same. The calculation of the decay rate gives, 

^ f'r. TTl^ 
J B s B s 



r NP (B s ^i+r) = gUgs+g'p" 1 ' 1 



a 



(29) 



The Branching ratio is given by, 

B(B S -> 1+1 



fkgHg's + g'p 



'2\ 



0.17 



(30) 



(m b + m s ) 2 

To estimate the value of g P {g's + g'p), we again consider the related decay B — » K*l + l~. Its 
matrix element, due to Lsp is given by, 



M(5 -> iTZ+r 



a: 



.2 

'W, 



# P (iT |s7 5 6| 5) g' s u{pi)v(pi) + g' P u(pi)~f 5 v(pi 



y/2 \ 47rs 

as (IT* |s6| 5) = 0. The pseudoscalar hadronic matrix element is given by [l 

/ 2rriK 



(31) 



(K* \s l5 b\ B) = -i 



\m b - m s 



(32) 



The q 2 dependence of the formfactor is described by a pole fit, 



A (q 2 



A)(0) 
;i - q 2 /m 2 B )' 



(33) 
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The full calculation gives, 



T NP (B -> K*l + l 



( 2rriK* 



2567T 3 J yAirsyy / \rrib — ra. 



where, 



Isp 



dz 



:i-z) 



i + 



m 



K* 



m, 



4m 2 ,, 
ml 



'SP 



(34) 



(35) 



The limits of integration for the dilepton invariant mass (z = q 2 /m 2 B ) are once again given 
by z min = (2mi/m,B) 2 and z max = (1 — m^./ms) 2 . Now we assume that the maximum value 
of this decay rate is the measured experimental value. Thus from Eq. (J3*4~j) . we get 



2 / '2 i '2 



'2N 



(m 6 - m s ) z B Exp {B - > K*Z+r' 
2.16L4 (0)] 2 



x 10 3 . 



(36) 



Taking the value of A (0) to be 0.47llg;^ (U|, we get 



^p(% 2 + ^p)=4.02™x10 



'2 



+2.41 



Substituting the value of g P (gs + fl'p) m Eq. (f30|) we get 



(37) 



B(B S -> /+r N 



2.10+J;| x 1(T 5 . 



(38) 



The upper bound on 5(-B s — > /i + M~) from the above equation is much higher than the 
present experimental upper bound Thus we see that the measured values of B(B — ► 
(K, K*)l + l~) do not provide any useful constraint on L 5 p contribution to B{B X — » 
The significance of this result is that if a future experiment, such as LHC-b [l6j observes 
B(B — > /! + /i~) > 10~ 8 , one can confidently assert that the new physics giving rise to 
this large a branching ratio must necessarily be of scalar /psuedo scalar type. Comparing the 
expression in Eq. (|30|) to the experimental upper bound in we obtain the bound 



9p(9s+9p)<W 



'2\ 



(39) 



Conclusions : We considered the most general effective Lagrangian for the flavour chang- 
ing neutral process b — > sl + l~, arising due to new physics. We showed that the present 
experimental values of B(B — > (K, K*)l + l~) set strong bounds on B(B S — > l + l~) if the 
effective Lagrangian is product of vectors/axial- vectors. Given that the above semi-leptonic 
decay rates of B-mesons are comparable to their SM predicted values, we showed that the 
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rate for purely leptonic decays of B s can't be much above the their SM predicted values. We 
have also derived a 3a upper bound on B(B S — > < 5 x 1CT 9 arising from Z'-mediated 

flavour changing neutral currents. If the effective Lagrangian for b — > sl + l~~ is product of 
scalars/psuedoscalars then present experimental values of B(B — > (K, K*)l + l~) do not lead 
any useful bound on B(B S — ► This leads us to the very important conclusion that, if 

a future experiment observes B s — > l + l~ with a branching ratio greater than 10~ 8 , then the 
new physics responsible for this decay must of be scalar/psuedoscalar type. 
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